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Abstract 

The JTZ model [C. Jung, T. Tel and E. Ziemniak, Chaos 3, (1993) 
555] , as a theoretical model of a plane wake behind a circular cylinder 
in a narrow channel at a moderate Reynolds number, has previously 
been employed to analyze phenomena of chaotic scattering. It is ex- 
tended here to describe an open plane wake without the confined nar- 
row channel by incorporating a double row of shedding vortices into 
the intermediate and far wake. The extended JTZ model is found in 
qualitative agreement with both direct numerical simulations and ex- 
perimental results in describing streamlines and vorticity contours. To 
further validate its applications to particle transport processes, the in- 
teraction between small spherical particles and vortices in an extended 
JTZ model flow is studied. It is shown that the particle size has signif- 
icant influences on the features of particle trajectories, which have two 
characteristic patterns: one is rotating around the vortex centers and 
the other accumulating in the exterior of vortices. Numerical results 
based on the extended JTZ model are found in qualitative agreement 
with experimental ones in the normal range of particle sizes. 

PACS numer(s): 05.45.-a, 47.32.-y 

Keywords: Plane wake, von Karman Vortex street, JTZ model. Particle 
dynamics 
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The plane wake behind a circular cylinder is one of the most fun- 
damental phenomena in fluid mechanics, involving periodic vortex 
shedding from the cylinder, which is known as von Karman vortex 
street. Recently, there emerge vast interests in the chaotic ad- 
vection of pairticles in the waike flow, in pcLrticuleir in its diversified 
applications to such processes as chemical and biological ones. The 
JTZ model was introduced to describe the plane wake in a narrow 
channel at a moderate Reynolds number and used to investigate 
chaotic advection of petrticles near the circuletr cylinder. In this 
paper we extend the JTZ model to describe an open plane wake 
without the confined narrow channel by adding a double row of 
shedding vortices in the intermediate and far wake, where the es- 
sential change is that Icirge damping in the original JTZ model is 
now replaced by little damping in the extended one. The results of 
extended JTZ model agree qualitatively with both direct numeri- 
cal simulations and experimental investigations on streamlines and 
vorticity contours. By using the extended JTZ model, particle 
transport processes £tre also simulated on different petrticle sizes. 
It is shown that the qualitative features of peirticle trajectories in 
the experimental investigation can be predicted by the numeri- 
cal simulation, justifying the application of extended JTZ model 
to quick and convenient estimations of the qualitative feature of 
peirticle transports in open plane weike flows. 
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1 Introduction 



The motion of particles in a non-uniform flow has received great attention 
due to its potential applications to both natural and engineering systems [1], 
as well as its dominant role played in the transport processes of particulate 
and multi-phase systems [2], such as those found in environmental engineer- 
ing, geophysical sciences, microfluids and combustion. In the transport pro- 
cesses, particles are not just passively carried by the background flow, but 
also have dynamics of their own. The motion of particles exhibits abundant 
characteristics even if the background flows are very simple. Particles tend 
to concentrate asymptotically along periodic, quasi-periodic or chaotic tra- 
jectories for such flows as steady and unsteady cellular flows[3,4], periodic 
Stuart vortex flows[5], the von Karman vortex street flows[6,7] and plane 
wake flows confined in a narrow channel (the JTZ model) [8,9]. 

The classical von Karman vortex street is a well known pattern [10, 11], 
and the viscous plane wake behind a circular cylinder, as an example of com- 
plex flow containing vortices and shear layers, has been extensively studied. 
For exhibiting the phenomenon of chaotic scattering, a theoretical model 
(the JTZ model) was proposed in [12] to describe plane wakes in narrow 
channels at a moderate Reynolds number (Re=250). It fits well the wake 
flow field given by direct numerical calculations and its application is much 
more convenient in getting the background wake flow field than solving the 
Navier-Stokes equations. The chaotic advection of particles near a circular 
cylinder has been investigated by using the JTZ model[13-17], however, rele- 
vant experiments of particle dispersion were conducted in open plane wakes 
consisting of near, intermediate and far regions[18,19]. And thus, we'll ex- 
tend in this paper the JTZ model to an open plane wake without the confined 
narrow channel. 
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2 Direct numerical simulation of a plane wake 



The incompressible plane wake is governed by the two-dimensional Navier- 
Stokes equations, as the following non-dimensionalized ones 

V«u = 0, 

(1) 

t = -Vp + iV(u) + ^L(u), 
where u is the fluid velocity and p is the fluid pressure divided by the fluid 
density. The uniform flow velocity Uoo and the diameter of circular cylinder 
D are taken as the characteristic variables of the system, and Re(= UooD/uf, 
I// is the fluids kinematic viscosity) is Reynolds number of the plane wake. 
In Eq. (1), the linear diffusion and nonlinear advection terms are described 

by 

^(u) = V^u, 

N{u) = -i[u • Vu + V • uu]. 
To solve the Navier-Stokes equations, we incorporate a high-order splitting 
algorithm into the spectral element method, please refer to [201 EI] for detail. 
Integrating the Navier-Stokes equations under the boundary conditions given 
in [21], we obtain velocity field u and vorticity field Uz at Re=250. The 
drag coefficient Cd and Strouhal number Sr are determined as 1.5 and 0.21, 
respectively, very close to the calculation results of [22| [23]. 

The streamlines and vorticity contours are displayed in Fig. 1, where a 
wave structure and the von Karman vortex street are shown in the down- 
stream wake with x restricted to the range of x < 20 to have a clear and 
regular vortex structure. Since the wake is a periodic process with a period 
of Tc = 1/Sr = 4.8, only figures at time phases and Tc/4 in the period Tc 
are plotted in Fig.l(a)(b) and Fig.l(c)(d), respectively. The corresponding 
figures at time phases Tc/2 and STc/A can be obtained by taking mirror im- 
ages of the plots of Fig.l(a)(b) and Fig.l(c)(d) about the x axis, respectively. 
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It is observed that only two vortices can be discerned in the streamhne plot 
of near wake. The regular von Karman vortex street appears in the inter- 
mediate and far wake as shown in the vorticity contour plots, where the 
semimajor axes of the elliptical vortices are in the vertical direction. Close 
to the cylinder, the distance between the vortices is about one quarter of 
those further away. Due to viscous decay, the vorticity decreases gradually 
in the wake downstream. The maximum absolute value of vorticity is very 
large in the boundary layer, and drops sharply in the near wake, and then, 
the maximum decreases slowly in the intermediate and far wake, having two 
constants in the range of 2 < x < 10 and 10 < x < 20, respectively. The 
maximum in the near wake is about 1.7 times that in the intermediate wake 
of 2 < X < 10, which is about 1.5 times that in the far wake of 10 < x < 20. 
In other words, the vortex structure in the wake downstream can be cap- 
tured in a snapshot and described approximately as follows. Within the 
wake region x < 20, there occur two transition zones with varying vorticity 
and two constant vorticity zones. The first transition zone is the near wake 
containing two vortices, and adjacent to which is the first constant vorticity 
zone, i.e. the intermediate wake with four vortices. And further downstream 
there occurs the second transition zone containing two vortices, and finally 
comes the second constant vorticity zone-the far wake with two vortices. An 
experimental investigation on open plane wakes at Re=250 was given in [21] 
with time phase around Tc/2, the streamline of which are in close agreement 
with the above-given simulation results. 
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3 An extended JTZ model 



The classical von Karman vortex street is the most simple model for a plane 
wake behind a circular cylinder, however, its computational zone doesn't 
contain the circular cylinder. The JTZ model for plane wakes docs cover the 
circular cylinder, and thus the flow field satisfies the no-slip boundary condi- 
tions at the cylinder surface. In what follows, we give a brief presentation of 
the JTZ model, where a plane wake behind a circular cylinder is considered 
within a narrow channel at a moderate Reynolds number (Re=250) with 
no-slip boundary conditions on the cylinder surface. By taking the cylinder 
radius Rq and the vortex shedding period Tc as the characteristic length and 
time of the flow, the dimensionless model stream function[12] is written as 



satisfies the no-slip boundary conditions at the cylinder surface. The second 
factor in Eq. (3) is 

g{x: y, t) = -whi{t)gi{x, y, t) + wh2{t)g2{x, y, t) + uoys{x, y). (5) 

The first two terms model the periodic detachment of the vortices, in which 
w represents average strength of the time-dependent vortices and 



^(a;,y,t) = f{x,y)g{x,y,t), 



(3) 



where the first factor 





hi{t) = |sin7rt|, 



(6) 



The factors 



giix,y,t) = exp-^o[(^-^iW)'+"'(2'-2'o)'l, 
g2{x,y,t) = exp-A[(^-^2{t))2+a2(2/+2/o)2] 



(7) 
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are the Gaussian forms with dimensionless vortex size /Jg with its center 
located at [xi{t),yQ] and [x2{t),—yo] in the wake. The vortices move down- 
stream at a constant velocity 

Xi(t) = 1 + Lomod()f:, 1), 
X2(t) = Xi(t - i). 

yo is the distance of the vortex centers from the x-axis, Lq is the dimensionless 
distance a vortex traverses during its lifetime and uq is the dimensionless 
background velocity. The last term in Eq. (5) arises from the background 
flow, and the screening factor 

s{x, y) = l- exp[-(a; - 1)V«' - y^] (9) 

ensures that the effect of the background flow velocity uq is reduced in the 
wake. In the numerical simulation, parameters are chosen as a = 1, a = 2, 
f3o = 0.35, Lq = 2, yo = 0.3, Uq = 14 and w = 24[15J. The streamlines and 
vorticity contours given by the JTZ model at time phases and 1/4 in a 
period are shown in Fig. 2, where the size is rescaled by 1/2 to set the length 
unit to D. 

Both the streamlines and vorticity plots show that the JTZ model gives 
only two vortices near the circular cylinder, which emerge and diminish in 
a period. So the JTZ model fits well the flow structures of plane wakes 
in a narrow channel. However, for open plane wakes without the confined 
narrow channel, the emerged vortices won't die out so quickly, and thus it 
can represent the flow structures in the near wake, rather than the structures 
of staggered double row vortices in the intermediate and far wake. To extend 
the JTZ model to the open plane wake, we will preserve the JTZ model for 
the near wake and add a double row of staggered vortices for the intermediate 
and far wake. 
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In the extended stream function ^(x, y, t) — f{x, y)g{x, y, t), the bound- 
ary function f{x,y) is kept, but g{x,y) is replaced by 

n 

g{x, y, t) = 70 Y.{-iyhi{t)gi{x, y, t) + uoys{x, y), (10) 

i=l 

where 

uo = UMR^, (11) 
gi{x, y, t) = exp-^o{[^-^*W]'+"?[2'-(-i)'"'2'ol'>, i = 1, 2, • • • , n. 

Here, n — 2 vortices are added in the intermediate and far wake flow field, 
and thus the extended model is more reasonable in representing the realistic 
structure of open plane wakes. F and C/qo are respectively the magnitude of 
vortex and free-stream velocity, and 70 and uq are their non-dimensionalized 
quantities based on Rq and Tg. The reference velocity is Rq/Tc instead of f/00, 
and the dimensionless inflow velocity is (ho, 0) instead of (1, 0). In Eq. (10), 
w of the original JTZ model is replaced by 70. The strength of vortices is 
defined as time-dependent in view of the viscous decay mentioned in section 
2, where n{= 10) vortices are observed in the range — 2 < x < 20 
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if mod{t, 1) < 0.5; 



(12) 



hi — |sm7ri|, 
h2 = l-0.8mod(t, 1), 
h^<i<& = 0.6, 
/i7 = 0.6 - 0.2mod(^, 1), 
hs = 0.5 -0.2mod(i, 1), 

^9<i<10 — 0.4, 

hi = 1 -0.8mod(f + 0.5, 1), 
h2 = |sm7r(t + 0.5)1 = |cos7rt|, 

^3<i<6 — 0.6, 

hy = 0.5 -0.2mod(i + 0.5,1), 
hs = 0.6 - 0.2mod(i + 0.5, 1), 
^9<j<io = 0.4, 

In the first half of the shedding period, the time dependent strength of a 
shedding vortex starts from and reaches 1. In the second half, it decreases 
from 1 to 0.6. In the later periods, it is changing or keeps unchanged depend- 
ing on its downstream distance. The vortex centers are located at positions 

[xi{t), ±yo], and move downstream with two different velocities Lo/2 and 2Lo 
in the first shedding period and the later periods, respectively. 



if mod(t, 1) > 0.5. 



xi{t) = l + ^mod(t, 1), 
X2{t) = l + ^mod(t + 0.5, i; 



Xi{t) 



Lo 



+ (i-3)Lo + 2Lomod(i,l), 



if i = 2A;-l(A;>2), 



1 + ^ + (i - 4)Lo + 2Lomod(i + 0.5, 1), if i = 2k{k > 2), 



(13) 



where Lq/A and Lq are the streamwise distance between two neighboring 
vortices within and outside the near wake region, respectively. Moreover, 
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s{x,y) in the last term of Eq. (10) is written as follows 

six,y) = 1 - exp[-(x - If/al - y% (14) 

In order to clarify essential features of the extended JTZ model, we display 
in Table 1 the time-variation of the strength and positions of four vortices 
near the cylinder in a period. In general, the up/down vortices in the double 
row of staggered vortices are denoted by odd/even. Each pair of vortices 
exists in a periodic streamwise interval [xi(=i^3^,,,)(t) , Xi+2(t)) , which is defined 
at time phase 0. For example, at time phase 0/|, vortex i = 1 and vortex 
i = 2 live in the first periodic streamwise interval between 1 and 1 + Lq/2, 
vortex i — 3 and vortex i — 4 live in the second one between 1 + Lo/2 and 
1 + 5Lo/2, and so on. After one half of the period, i.e., at time phase |/0, a 
new down/up vortex is shed from the circular cylinder and denoted by 2/1. 
At the same time, the down/up vortices in the double row of vortices move 
away from the original periodic streamwise positions and are re-located at the 
next periodic streamwise positions. Therefore, the notation i for the original 
down/up row of vortices at time phase 0/| will be rewritten as i + 2. For 
example, ait = 1/2, vortex i = 2 and vortex i = 4 are located at x = l-|-Lo/2 
and l + 5Lo/2, respectively, i = 4 and i = 6 will be used to denote the vortex 
i = 2atx = l-|-Lo/2 and vortex i = 4atx = l-|- 5Lo/2, respectively. At 
the same time, a new vortex is shed at (1, —yo). It will be denoted by i — 2. 
At t — 1, vortex i — 1 and vortex i — 3 are located at x — 1 + Lo/2 and 
1 -|- 5Lo/2, respectively, i = 3 and i = 5 will be used to denote the vortices 
i = lata; = l-|- Lo/2 and i = 3atx = l + 5Lo/2, respectively. At the same 
time, a new vortex is shed at (1, yo)- It will be denoted by i 1. The vortex 
shedding process is periodic with period t — 1. 

Parameters n — 10, a — 1, ao — 2, Po — 0.35 and yo = 0.3 are used in the 
following simulations and parameter ctj is defined as follows 
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2 - 1.2mod(i,l), if i = 1; 
a,; = < 2 - 1.2mod(i + 0.5, 1), if i = 2; (15) 
0.8, if i > 3, 

which describes the continuous change of the horizontal axes of eUiptical 
vortices from the semimajor in the near wake to the semiminor in the in- 
termediate and far wake. Prom the definition of Sr{— 2i?o/C/oo^c), the non- 
dimensional velocity Uq is written as Uq = 2 /Sr. For Re = 250, Sr equals 0.21 
as given in [22]. Based on the velocity and vorticity fields given by the direct 
numerical simulation described in section 2, we can determine approximately 
7o = 2.2uq and Lq = 4.4. 

The streamlines and vorticity contours at time phases and 1/4 in a 
period are shown in Fig. 3, which is better than the original JTZ model in 
capturing the global distribution of vortices in the open plane wake field given 
by the direct numerical simulation. In particular, the vorticity distribution 
in the intermediate and far wake of the extended JTZ model is found in 
quantitative agreement with that of direct numerical simulation, except for 
the region near the cylinder surface where there exists some difference. In 
the direct numerical simulation, the vorticity is confined to a thin boundary 
layer, whereas in the extended JTZ model, it is extended to a finite domain 
adjacent to the cylinder surface. Anyway, the extended JTZ model can 
be conveniently used as a simple plane wake model to investigate particle 
dynamics in plane wake fiows. 



4 Particle dynamics in the extended JTZ model 

For the motion of a small spherical particle in fluid flow, it is well-known 
that the flow fleld around the particle can be approximated by a Stokes flow 
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provided the particle diameter is small enough compared to the characteristic 
length of the fluid motion. Therefore, the particle-particle interactions as well 
as the effects of particles on the flow are negligible. Under these assumptions, 
the momentum equation of the motion of the small spherical particle [2 5] is 
written as follows 



(16) 

where V is velocity of the particle, and p is density, g is gravitational acceler- 
ation, and subscripts / and p refer to the fluid and particle, respectively The 
derivatives D/Dt and d/dt are used to denote the time derivative following 
a fluid element and the moving sphere, respectively. Introducing the dimen- 
sionless quantities 6 — Pp/pf, e — 1/(0.5 + 5), x* = x/i?o, t* — t/{Ro/U^), 
u* = u/C/oo, V* = V/C/oo and g* = g/g, based on i?o and U^o, we can 
non-dimensionlize Eq.(16) as follows 

dV ^ g-l-f )g I 3e Du I ^/ yN 

dt Fr^ & ^ 2 Dt ^ Stv" ^/ ^-^rj-^ 

where Fr = Uoo/ V qRq and St = UooT/Rq {T is the particle viscous relax- 
ation time, dp/lSeuf) are the Proude number and the Stokes number, respec- 
tively. The density and size of particles in Eq.(16) are taken from [18] as 
Pp = 2.4 X lO^kg/m^ and dp = 0(10~^)m, respectively. Since air is chosen as 
the fluid medium in the flow, the fluid properties in Eq.(16) are described by 
Pf = 1.225kg /m^ and Uf — 1.45 x 10~^m^/s[26]. To emphasize the effects of 
viscous force or St in Eq. (17) on particle transport process, the background 
flow field parameters are chosen as U^o — 0.1m/ s and Rq — 1.81 x 10~^m. 
Since S is fixed and is of the order of O(IO^), parameter e appears to be of 
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the order of 0(10~^). The particle diameters dp is very small [dp — 0(10~^)m], 
so the parameters 1/St and appear to be of the order of O(IO^) and 

0(10°), respectively. In this case, we can neglect the following smaller order 
terms: the gravity term, the stress tensor term, the Basset history term and 
reduce the equation of motion (17) to 

In section 3, the physical quantities are non-dimensionalized using and 
Rq, which are different from the above ones. The function u{x,y,t) in Eq. 
(18) is now replaced by 

u(a;, y, t) = u{x, y, i)/uo, 

where i and u are the time and velocity field for the extended JTZ model, 
respectively. In the definition of Sr, the characteristic length is 2Rq, rather 
than i?o, which is the characteristic length of the model. Since i x Tc — 
t X Rq/Uoo, so i ^ t X Ro/{UooTc) ^ t x Sr/2. In what follows, we will 
analyze essential features of particle trajectories in the extended JTZ model 
fiow. 

Equation (19) is integrated with time step At = 0.001 — 0.01 by using a 
fourth-order Runge-Kutta algorithm. Particles are released along a semicircle 
(of radius 1.2) with time interval 0.05. In each time interval, the released 
particles total to 100. The initial velocities of particles are taken as the local 
fiow ones. We choose t = as the time for the initial release of particles, and 
give a snapshot of particle distribution in the wake at ^ = 3. 

The particle diameter dp is taken in the range of 8 x 10~^ — 8 x 10~^m, 
which corresponds to St — 0.003 — 0.27. In the numerical simulation of 
dp = 8 X 10-^m, 1 X 10-^m, 2 x lO^^m, 3 x lO^^m, 5 x lO^^m and 8 x lO^^m, 
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the particles rotate around neighboring vortices as well as move downstream 
under the effects of viscous force (u — V) / St. In interacting with vortices, 
they exhibit some characteristic distributions. According to the rotation di- 
rections of vortices, the particle clusters form foliations near the vortices, and 
show two typical patterns depending on the particle sizes. Firstly, when the 
particle diameter is small enough (in the range of 8 x 10~^ — 2 x 10~^m), the 
foliations rotate around the vortex centers with an example shown in Fig. 
4(a) for dp = 1 X 10~^m [St = 0.004). The large viscous force provides cen- 
tripetal force to particles, and thus the foliations are folded by the rotation of 
the vortices and stretched by the training field. When the particle diameter 
is in the range of 2 x 10~^ — 8 x 10~^m, the foliations are accumulated in the 
exterior of vortices with an example shown in Fig. 4(b) for dp = 5 x lO^^m 
{St — 0.1). Since the viscous force decreases as dp increases, the larger the 
dp, the smaller the centripetal force provided by the viscous force. Hence, 
the foliations are moved away from the vortex centers and accumulate in 
the exterior of vortices. In the experiments of [18,19], the particle dispersion 
in a plane wake is investigated in detail. As the particle size increases, the 
organized pattern of particle trajectories changes from filling in the vortex 
structures to accumulating along the exterior of vortices. In the direct nu- 
merical simulation of [27] , the particle distributions for different particle sizes 
are similar to the above-mentioned experimental observations, and these be- 
haviors agree qualitatively with our numerical results based on the extended 
JTZ model. 
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5 Conclusion 



In summary, we generalize and extende the JTZ model, which is a theoret- 
ical model for a plane wake behind a circular cylinder in a narrow channel 
at a moderate Reynolds number, to describe an open plane wake without 
the confined narrow channel by incorporating a global distribution of shed- 
ding vortices. The extended JTZ model is found in quahtative agreement 
with both direct numerical simulations and experimental studies in respect 
of streamlines and vorticity contours. The interaction of small spherical par- 
ticles with vortices in the extended JTZ model flow is investigated to further 
validate the moderling in studying particle transport processes. It is found 
that the particle size has significant infiuences on the features of particle tra- 
jectories, and the particles exhibit two typical patterns as the particle size 
increases. When the particle size is small enough, particle clusters would 
rotate around the vortex centers, whereas they would accumulate in the ex- 
terior of vortices for larger particle size. The numerical results are found in 
qualitatively agreement with the experimental ones in the normal range of 
particle sizes. 

Acknowledgments The author thanks ICTS and IMECH research com- 
puting facilities for assisting in the computation. 
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Table 1. Time- variation of the strength and positions of four vortices near 
the cyhnder in a period. 



t 





1/4 


1/2 


3/4 


hi 





V2/2 


1 


0.8 


h2 


1 


0.8 





V2/2 


hs 


0.6 


0.6 


0.6 


0.6 


hi 


0.6 


0.6 


0.6 


0.6 


Xi 


1 


l+Lo/8 




I+3L0/8 


X2 


l+Lo/4 


I+3L0/8 


1 


l+Lo/8 


X3 


l+Lo/2 


1+Lo 


I+3L0/2 


I+2L0 


X4 


I+3L0/2 


l+2Lo 


l+Lo/2 


1+Lo 
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FIGURE CAPTION 

Fig. 1. Streamlines and vorticity contours for the plane wake flow at Re — 
250 obtained from direct numerical simulation. The values of vorticity ojd,u^, 
which are non-dimensionalized by D and U^, are shown in the vorticity 
contour plots. The time phases are chosen as in (a)— 2 < x < 10; (b) 
10 < X < 20 and T^/A in (c)-2 < x < 10; (d)10 < x < 20 in a vortex 
shedding period Tc. 

Fig. 2. Streamlines and vorticity contours for the JTZ model flow at time 
phases (a) and (b) 1/4 in a period. 

Fig. 3. Streamlines and vorticity contours for the extended JTZ model 
flow at time phases in (a) -2 < x < 10; (b) 10 < x < 20 and 1/4 in (c)-2 < 
X < 10; (d)10 < a; < 20 in a period. The index i of vortices in the extended 
JTZ model (10) is marked at the upper and bottom zones of the vortex street 
central region. The non-dimensional values ujd,Uoo the vorticity contour 
plots are calculated from the modelling of vorticity ijJrq^To scaled by Rq and 
Tc by referring to the relation ojd,u^ — Stujjiq^Tc = 0.21ci;rq ^c- 

Fig. 4. Instantaneous particle distribution and vorticity contour in the 
extended JTZ model flow at time phase in a period with different particle 
diameters (a) dp — lx 10~^m and (b) dp — 5x 10~^m. The release of particles 
starts at t — and is completed at t — 3 when the snapshot is taken. 
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